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High-energy photons from pair annihilation of dark matter particles contribute to the cosmic 
gamma-ray background (CGB) observed in a wide energy range. Since dark matter particles are 
weakly interacting, annihilation can happen only in high density regions such as dark matter halos. 

The precise shape of the energy spectrum of CGB depends on the nature of dark matter particles — 
their mass and annihilation cross section, as well as the cosmological evolution of dark matter halos. 

In order to discriminate between the signals from dark matter annihilation and other astrophysical 
sources, however, the information from the energy spectrum of CGB may not be sufficient. We show 
that dark matter annihilation not only contributes to the mean CGB intensity, but also produces 
a characteristic anisotropy, which provides a powerful tool for testing the origins of the observed 
CGB. We develop the formalism based on a halo model approach to analytically calculate the three- 
dimensional power spectrum of dark matter clumping, which determines the power spectrum of 
annihilation signals. We show that the expected sensitivity of future gamma-ray detectors such as 
the Gamma Ray Large Area Space Telescope (GLAST) should allow us to measure the angular 
power spectrum of CGB anisotropy, if dark matter particles are supersymmetric neutralinos and 
they account for most of the observed mean intensity of CGB in GeV region. On the other hand, if 
dark matter has a relatively small mass, on the order of 20 MeV, and accounts for most of the CGB 
in MeV region, then the future Advanced Compton Telescope (ACT) should be able to measure the 
angular power spectrum in MeV region. As the intensity of photons from annihilation is proportional 
to the density squared, we show that the predicted shape of the angular power spectrum of gamma 
rays from dark matter annihilation is different from that due to other astrophysical sources such as 
blazars and supernovae, whose intensity is linearly proportional to density. Therefore, the angular 
power spectrum of the CGB provides a “smoking-gun” signature of gamma rays from dark matter 
annihilation. While the mean GGB intensity expected from dark matter halos with smooth density 
profiles is smaller than observed, the dark matter substructure within halos may provide the origin 
of additional “boost” factors for the annihilation signal. Our formalism can be used for any other 
radiation processes that involve collision of particles. 

PACS numbers: 95.35.-|-d; 95.85.Pw; 98.70.Vc 


I. INTRODUCTION 

The energy density of the universe is dominated by 
invisible components: dark matter and dark energy, both 
of which are of unknown origin. Remarkable progress 
in both the theoretical and observational studies of the 
cosmic microwave background (CMB) anisotropy. Type 
la supernovae, and large-scale structure has allowed us 
to precisely determine what fraction of the total energy 
these dark components convey, = 0.23 and UIa = 0.73 
Q, where the subscripts x A denote dark matter 
particles and dark energy, respectively. 

While the nature of dark energy is a complete mystery, 
there are several candidates for dark matter particles. Of 
which, the most popular candidate is the stable super- 
symmetric neutralino with mass on the order of GeV to 
TeV, which can explain the observed dark matter mass 
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density today. Although very weak, it is expected that 
dark matter particles may interact with the usual mat¬ 
ter via scattering, interact with themselves, and/or an¬ 
nihilate into gamma rays, positrons, and neutrinos. The 
direct detection of particle dark matter is, therefore, un¬ 
der intensive and extensive efforts of both physicists and 
astronomers iii 

High-energy photons from annihilation of dark mat¬ 
ter particles provide indirect means to probe the proper¬ 
ties of dark matter. Annihilation signatures, especially 
gamma rays, have been searched for in regions where the 
dark matter density is expected to be high, as annihila¬ 
tion rate is proportional to the density squared, p^. An 
obvious site is the central region of our Galaxy. Strong 
gamma-ray emission has been detected toward the Galac¬ 
tic center over a wide energy range, and the nature and 
origin of this emission have been investigated by many re¬ 
searchers 0 IE 0 B Q • Although both the 

spectrum and angular distribution of gamma rays from 
the Galactic center are well observed, it is still quite dif¬ 
ficult to distinguish a dark matter component from the 
other possibilities such as emission from “ordinary” as¬ 
trophysical objects; even future gamma-ray detectors and 
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telescopes do not have sufficient resolution to remove all 
the other sources from the Galactic central image. 

Another possibility is the extragalactic background 
light, the cosmic gamma-ray background (CGB), which 
has been measured in a wide energy range [ij, ll^ lift 
in in It has been speculated that some fraction of 
the CGB may originate from annihilation of dark mat¬ 
ter particles in halos distributed over cosmological dis¬ 
tances in m m m El H El Es m- The dark 

matter contribution to the CGB depends on the nature 
of dark matter particles, such as their mass and anni¬ 
hilation cross section, as well as the cosmological evo¬ 
lution of dark matter halos. While it is likely that the 
dominant part of the CGB spectrum in the GeV region, 
detected by the Energetic Gamma Ray Experiment Tele¬ 
scope (EGRET) mmij comes from unresolved blazars, 
i.e., beamed population of Active Galactic Nuclei (AGN), 
annihilating dark matter might give large contribution at 
some specific energy range.^ Possible candidates for this 
case include the supersymmetric neutralino, as well as 
the Kaluza-Klein dark matter medicted by theories of 
universal extra-dimension (e.g., y) and heavy relic neu¬ 
trinos 1^ Is^ ■ On the other hand, the origin of the CGB 
in MeV region is much less understood than in GeV re¬ 
gion. The soft gamma-ray spectrum in 1-20 MeV cannot 
fully be attributed to either AGN or Type la supernovae 
or a combination of the two |3lLl32| : thus, annihilation of 
dark matter particles is one of the most viable explana¬ 
tions for the CGB in MeV regi on, p rovided that the dark 
matter mass is around 20 MeV |2^|^. Such a light dark 
matter particle was originally introduced to explain the 
origin of the 511 keV emission line from the Galactic cen¬ 
ter 1^ , detected by the International Gamma-Ray 
Astrophysics Laboratory (INTEGRAL) |^ |^ , which is 
otherwise difficult to explain. Therefore, the MeV dark 
matter is an attractive candidate that satisfies the obser¬ 
vational constraints from both the Galactic center and 
the CGB,^ while its particle physics motivation is less 
clear than for the neutralinos. 

Dark matter annihilation may be a viable explanation 
for the CGB, but how do we know for sure that the CGB 
does come from annihilation? What would be a smoking- 


^ As has been known and is shown below, with the canonical 
choice of relevant parameters, dark matter annihilation cannot 
give enough contribution to the mean CGB intensity; one still 
needs large boost by some mechanism. In addition, its contribu¬ 
tion is strongly cons trai ned by gamma-ray observations toward 
the Galactic center 1^ . However, one might avoid these con¬ 
straints by invoking gamma-ray em ission from sub-halos having 
M ~ 10“®Mq within host halos |2^, or fro m d ensity cusps form¬ 
ing around intermediate-mass black holes |2^ 

^ After this paper has been submitted. Ref. [23 has appeared on 
the preprint server. They claim that the mass of MeV dark mat¬ 
ter should be less than 3 MeV from a more accurate treatment 
of relativistic positrons annihilating with electrons in the inter¬ 
stellar medium. While we use = 20 MeV throughout this 
paper, one can easily extend our calculations to arbitrary dark 
matter masses. 


gun signature for the annihilation signal? We argue that 
anisotropy of the CGB may provide a smoking-gun sig¬ 
nature. Although the CGB is isotropic at the leading 
order, anisotropy should also exist if the CGB originates 
from cosmological halos. The future gamma-ray detec¬ 
tors with an enhanced sensitivity and angular resolution, 
such as the Gamma Ray Large Area Space Telescope 
(GLAST) or the Advanced Compton Telescope (ACT), 
should be able to see such anisotropy. We believe that 
the CGB anisotropy is going to be the key to discrimi¬ 
nating between the dark matter annihilation signal and 
the other sources. 

In this paper, we develop the formalism to analyti¬ 
cally calculate the angular power spectrum of the CGB 
from dark matter annihilation. We calculate the angular 
power spectrum in GeV region (for supersymmetric neu¬ 
tralinos) as well as in MeV region (for MeV dark matter). 
We then discuss the detectability of CGB anisotropy in 
GeV region by GLAST and in MeV region by ACT, show¬ 
ing that the predicted anisotropy can be easily measured 
by 1-year operation of these experiments. The formal¬ 
ism given in this paper can also be used to evaluate the 
CGB anisotropy due to the other astrophysical objects 
with some modification. Note that the case of Type la 
supernovae has been discussed in Ref. Is^, being com¬ 
plementary to the present study. 

This paper is organized as follows. In Sec. QH we review 
the mean intensity of the CGB from dark matter annihi¬ 
lation. In Sec. ITTTI we develop the formalism to calculate 
the power spectrum of CGB anisotropy in the context of 
the halo models. Results of the angular power spectrum 
are then shown in Sec. CYI for both the neutralino and 
MeV dark matter. We compare the predictions with the 
expected sensitivities of the future gamma-ray detectors. 
In Sec. m we discuss energy dependence of the obtained 
anisotropy (Sec. WJh . other astrophysical sources which 
might contaminate the cosmological CGB (Sec. IVBll . and 
the effects of dark matter substructures (Sec. IWl . In 
Sec. EH we conclude the present paper with a brief sum¬ 
mary. 


II. COSMIC GAMMA-RAY BACKGROUND: 

ISOTROPIC COMPONENT 

Since the gamma-ray emissivity from dark matter an¬ 
nihilation is proportional to the density squared, the 
CGB intensity (the number of photons per unit area, 
time, solid angle, and energy range) toward the direction 
n can be generally expressed as 

I^{'h,Ej)= J dr S'^{r,n.r)W{[l + z]E^,r), (1) 

where E-y is the observed gamma-ray energy, z is the 
redshift, r is the comoving distance out to an object at 
z, IV is some function of gamma-ray energy and r that 
is given below, and 5 is the overdensity at hr compared 
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to the average dark matter mass density, given by 


The clumping factor, (d^), is then given by^ 


5=^. (2) 

Px 

Because the dominant contribution comes from the dark 
matter halos, (5 is always much larger than unity. Note 
that we label time by r (or redshift z used interchange¬ 
ably), and space by r = nr. 

We first derive the specific form of the function W. 
The general formula for the intensity is given by 


EjI.y{n,Ery) = 


_£_/■, P.y([l -f z]E^,z,nr) 
4 tt J ^ H{z){l + z)"^ 


X e 


— t{[1+z]E^,z) 


( 3 ) 


where P.y is the volume emissivity (energy of photons 
per unit volume, time, and energy range), H^{z) = 
HQ[flrn{^ + is the Hubble function in a flat uni¬ 

verse, and we assume the standard values for cosmologi¬ 
cal parameters, Hq = 100 h km s“^ Mpc“^ with h — 0.71, 
= 0.27, and JIa = 0.73. By the exponential factor, 
we incorporate the effect of the gamma-ray absorption 
due to pair creation with the diffuse extragalactic back¬ 
ground light in the infrared or optical bands how¬ 
ever, this effect is actually negligible for the cases we con¬ 
sider here where the gamma-ray energy is smaller than 1 
TeV. The annihilation rate per unit volume is given by 
{p^/m^)‘^{av)l2, where is the mass of dark matter 
particle and {av) is the annihilation cross section times 
relative velocity averaged with the weight of the velocity 
distribution. The volume emissivity is therefore given by 


P^{E^,z,nr)=E^^- 


Px{z,nr) 


rriy 


( 4 ) 


where dN^/dE^ is the gamma-ray spectrum per anni¬ 
hilation. Comparing Eq. m with Eqs. m and 0 , 
and recalling that dr = —cdz/H{z) and py^(z,hr) = 
p^{z)d{z,hr) = Q^pc{l + z)^S{z,hr), we obtain W as 


W{E^,z) 


Stt 


TOx ) 


2 


(1 + Z)3 


dEj 


( 5 ) 


By taking the ensemble average of Eq. o, we ob¬ 
tain the isotropic background component of the inten¬ 
sity, {Ij{Ej))-, it is equivalent to our evaluating {6‘^{z)), 
the mean clumping factor of dark matter density fluctua¬ 
tions, which is given by the integral of two-point correla¬ 
tion over the momentum. We use a halo model approach 
(see Ref. ^3 fo'' ^ review) to calculate a reduced part of 
the A^-point correlation function of density fluctuations. 





0,71 


X / dV w^(r|M, z). 



( 6 ) 


where u{r\M, z) = p(r|M, z)/M is the density profile di¬ 
vided by the halo mass M, dn/dM the halo mass function 
for which we use the expression given in Ref. This 
clumping factor (d^(z)) is also known as an intensity mul¬ 
tiplier (it is the same quantity as, e.g., A^(z) of Ref. pn| 
and f{z) of Ref. [^1. 

We should comment on the boundary of the mass in¬ 
tegral in Eq. ®. As the mass function falls exponen¬ 
tially toward large masses, the integral converges at high 
mass end; however, it actually diverges at low mass end, 
although only very weakly. This raises some questions 
regarding robustness of theoretical predictions for the 
mean intensity of CGB from dark matter annihilation. 
What determines the minimum mass, Mmin? In prin¬ 
ciple, one should integrate the mass function down to 
the free-streaming scale or the Jeans mass of dark mat¬ 
ter particles, below which there is no fluctuation in dark 
matter. For both neutralinos and MeV dark matter, the 
minimum mass is roug hly on the order of the mass of 
Earth [il liM li^ while it can vary by orders of 
magnitude depending on the precise value of the dark 
matter mass. (One should recalculate Mmin as a func¬ 
tion of the dark matter mass or any other parti¬ 
cle physics parameters that affect the strength of inter¬ 
actions.) When gamma rays emerge from annihilation 
without any other interactions of the products of annihi¬ 
lation with the materials within halos — i.e., dark mat¬ 
ter annihilating directly into photons, photons emerge as 
the final state radiation of annihilation (also known as 
the internal bremsstrahlung), or decay of the products of 
annihilation into gamma rays (e.g., decay of tt^ produced 
by annihilation of neutralinos), one should use the min¬ 
imum mass set by the free-streaming mass or the Jeans 
mass of dark matter particles. On the other hand, when 
gamma rays are produced predominantly by interactions 
between the products of annihilation and the materials 
within halos after annihilation, e.g., dark matter annihi¬ 
lation producing intermediate particles such as positrons 
which, in turn, annihilate with electrons in the halo to 
produce gamma rays, the minimum mass would be set 
by the Jeans mass of baryons, which is on the order of 
10^Mq or larger ^3. Therefore, which minimum mass 
to use depends on what radiative processes one assumes 
for the production of gamma rays per annihilation. 


^ The expression is identical to the 1-halo term of 

Eq. (2.15a) in Ref. M- The second (2-halo) term disap¬ 
pears, because we evaluate the quantity at one point and halos 
are assumed to be exclusive, having stiff boundary. 
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This issue is further complicated by the fact that not all 
dark matter halos would actually survive during the for¬ 
mation of large-scale structure. Halos above the baryonic 
Jeans mass would survive without any problems; how¬ 
ever, a signihcant fraction of halos with very small free- 
streaming masses, say 10 “®M©, would have been tidally 
disrupted shortly after their formation and might not be 
able to contribute to the CGB. The critical mass above 
which the tidal disruption due to hierarchical clustering 
becomes inefficient is difficult to calculate accurately at 
present [4^|^. (Note that we are not concerned about 
the tidal disruption of these micro halos within larger 
halos — this would be a subject of the effect of substruc¬ 
tures within halos, and we shall come back to this point 
in Sec. \VU Here, we are mostly concerned about the 
tidal disruption of the extragalactic micro halos that are 
outside of larger halos.) 

Generally, the dependence of GGB anisotropy on the 
minimum mass is expected to be not very strong, when 
the amplitude of anisotropy is divided by the mean in¬ 
tensity, < 5 / 7 /(Ty), the dependence approximately cancel¬ 
ing out in the numerator and denominator. To take into 
account the uncertainty in our understanding of the min¬ 
imum mass contributing to the GGB, we explore two 
cases, M„iin = 10 ®M© and 10 “®M©, i.e., 12 orders of 
magnitude variation in the minimum mass. 

The clumpiiig factor depends on z through the halo 
mass function Id.Jl and the concentration of density pro¬ 
files 1^1^ . To compute the linear power spectrum that 
is necessary for the halo mass function, we adopt the 
fitting formula in Ref. [bill . The concentration param¬ 
eter as a function of halo mass and redshift determines 
how steeply dark matter distributes in a halo; we use the 
result given in Ref. [^. Using these ingredients, we com¬ 
pute the clumping factor for a given density profile. For 
the Navarro-Frenk-White (NFW) profile (7 = !> 

where 7 is defined by p oc r~'^ for small radii), we obtain 
((52(0)) ~ 4 X 10^ for Mmin = 10®M© and 3 x 10® for 
Mjnin = 10“®M©. For steeper profiles such as the one 
with 7 = 1.5 suggested by the other simulations |^, the 
clumping factor increases by about an order of magni¬ 
tude. More recent simulations suggest that 7 depends on 
radii and could be anywhere between 1 and 1.5 in the 
inner region (see also Refs. 1^ 113 

analytical approaches). Again, while the mean intensity 
depends strongly on 7 (see, e.g.. Ref. [1^), anisotropy 
(diyided by the mean intensity) is expected to be less 
sensitiye to the precise yalue of 7. For definiteness we 
adopt 7=1 throughout the paper. Incidentally, it is 
practically useful to use the NFW prohle, as many quan¬ 
tities used in the formalism can be calculated analytically 
for this profile. 

One also needs to specify the particle physics parts 
inyolved in the function W — the dark matter mass, 
averaged cross section (av), and gamma-ray spectrum 
per annihilation. We assume that the dark matter par¬ 
ticle contributing to the GGB totally accounts for the 
observed dark matter mass density, = 0.23. In many 


cases, the annihilation cross section is closely related to 
as it determines the abundance of dark matter par¬ 
ticles having survived pair annihilation in the early uni¬ 
verse.^ For both the supersymmetric neutralino and MeV 
dark matter, the canonical value of the cross section is 
{av) = 3 X 10“26 cm® s“®, while a wide range of param¬ 
eter space is still allowed®; in addition, we assume it to 
be independent of the relative velocity. Once again, the 
value of the annihilation cross section actually does not 
affect the GGB anisotropy, when it is normalized by the 
mean intensity, SI^/{Ij). (The cancellation is exact when 
it is independent of velocity.) Thus, although a larger 
cross section is favored in order to make the dark matter 
contribution to the GGB mean intensity more signihcant, 
the prediction for the amplitude of anisotropy (divided 
by the mean intensity) is robust regardless of the precise 
value of the cross section. 

While the energy spectrum of gamma-ray emission per 
annihilation contains rich information about the nature 
of dark matter particles, we do not pay much attention 
to the detailed shape of the spectrum, but rather simply 
split the spectrum into two parts: MeV and GeV regions, 
the former representing MeV dark matter, and the latter 
representing neutralinos. This is a reasonable approxima¬ 
tion, as each contribution gives continuum gamma-ray 
emission that is roughly constant over relevant energy 
regions. For neutralinos we use a simple parameteriza¬ 
tion, ~ (0.73/TO^)e“'^-’^®-®'^/’"’(/[(U.y/m^)®-® -I- 

0.00014] ^ 3 . The Kaluza-Klein dark matter parti¬ 
cle might also contribute to the GGB in GeV re¬ 
gion; their contribution may be calculated using an 
appropriate gamma-ray energy spectrum per annihila¬ 
tion. For MeV dark matter we consider the internal 
bremsstrahlung (yy ^ e+e“ 7 ) as the source of contin¬ 
uum emission, which gives dN^/dE^ = a\\'n.{s'/rn^) — 
1][1 -I- (s'/s)2]/(7rU-y), where a = 1/137 is the fine struc¬ 
ture constant, s = and s' = Am^{m^ — E^) 

We shall ignore the 511 keV line emission, as its contri¬ 
bution must be subdominant compared with the AGN 
contribution 

In Fig. n we show the average GGB intensity cal¬ 
culated with the ingredients given above for the neu¬ 
tralino {rriy, = 100 GeV) and MeV dark matter (to^ = 20 
MeV). We also show the observational data obtained by 
the Imaging Gompton Telescope (GOMPTEL) and 


^ This would not be quite true when co-annihilations of dark mat¬ 
ter particles with other particles dominate (e.g., neutralinos co¬ 
annihilating with other supersymmetric particles). In such a case 
the annihilation cross section would have little to do with 
See Ref. for more details. 

^ For MeV dark matter, it has been argued that the canonical, 
velocity-independent cross sec tion would overproduce the Galac¬ 
tic 511 keV emission l33L l63l . However, uncertainty in our un¬ 
derstanding of the dark matter profile of our Gal axy stil l allows 
for the canonical cross section for > 20 MeV l26l l27ll . If, on 
the other hand, turns out to be less than 20 MeV, then a 
smaller cross section may be preferred. 
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E, [GeV] 


FIG. 1: Average intensity of the CGB from dark matter an¬ 
nihilation compared with the COMPTEL and EGRET data. 
The curve in GeV region is the predicted intensity for neu- 
tralinos with = 100 GeV, while that in MeV region is the 
prediction for MeV dark matter particles with = 20 MeV. 
As for the normalization, we have used {av) = 3 x 10“^® cm® 
s“^, and extra “boost” factors have been multiplied to match 
the observational data. Eor details of the boost factors, see 
the last paragraph of Sec. II. 


EGRET M The predicted spectra have been multiplied 
by a factor of 1.5 x 10® (2.4 x 10^) for the neutralino, and 
10 (1.6) for the MeV dark matter, respectively, assum¬ 
ing Mmin = 1 O®M 0 (IO^^Mq), to reasonably match the 
observational data. These boost factors may originate 
from dark matter substructures within halos or the cross 
section potentially larger than the value adopted here. 
Note that the 20 MeV dark matter with Mmin = 10“®Mq 
(which is not very far away from the actual free-streaming 
mass of 20 MeV dark matter particles) can already ac¬ 
count for most of the observed CGB without the need for 
a substantial boost due to substructures (the boost factor 
is only 1 . 6 ), which ag rees with a more detailed analysis 
performed in Ref. |27l| . 


III. COSMIC GAMMA-RAY BACKGROUND 
ANISOTROPY 

In this section, we develop the formalism to analyt¬ 
ically calculate the angular power spectrum of CGB 
anisotropy due to annihilation of dark matter particles. 
We note that this formalism is quite general, and with 
some modification, it can be applied to other situations 
including gamma rays from ordinary astrophysical ob¬ 


jects. 

A. General setup 

We expand the deviation of CGB intensity from its 
mean value into the spherical harmonics coefficients, aim, 

— Ty(?T.) — (-^ 7 ) ^ ^ aim^lmij^^- (7) 

Im 

Note that we have defined aim as a dimensionless quan¬ 
tity; thus, aim represents the amplitude of anisotropy 
divided by the mean intensity. Using the orthonormal 
relation of Yim{n), one obtains 

{Ij)aim = J dn 

where 

(9) 

and the energy index in the function W has been sup¬ 
pressed to simplify the notation. 

The goal of the present paper is to evaluate the angular 
power spectrum, C/ = (|a;mP), which is given by 

il7)"Ci= J ^ {W{[l + z]E^,r)}^Pf(^k=pTy 

( 10 ) 

The detailed derivation of this result is given in Ap¬ 
pendix El Here, Pf(k) is the three-dimensional (3D) 
power spectrum of /, defined by 

(/fc/fc') = (27r)®5(®)(/c -f k')Pf{k), (11) 

where represents the A-dimensional delta function, 
and fk is the Fourier transform of /(r). Note that k 
is a comoving wavenumber, as it is a Fourier variable 
corresponding to the comoving distance r. Given this 
relation m, we first focus on deriving the 3D power 
spectrum, Pf(k), and then turn to evaluating the angular 
power spectrum. 

B. Three-dimensional power spectrum 

1. Two-point contribution 

The power spectrum, Pf(k), is the Fourier transform 
of the two-point correlation function of / in real space, 

/n\ 

Ef \x — y) = {f{x)f{y)). As / is a quadratic function 
(‘ 2 ') 

of 6, Ef ' is given by the two- and four-point correlation 
functions of <5, and ^ 

^f\x-y) = S!^‘^\x,x,y,y) + 2S}-'^\x-yf, {12) 
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FIG. 2: Dimensionless power spectrum oi f — 5^ — (S^), 
A^ 2 (^) [Eq. 115ll 1. from the two-point correlation term 
[Eq. 114ll 1. evaluated at z = 0 for Mmin = 10® Mq. 


where we have suppressed the label r (or equivalently z) 
to simply the notation. This motivates our decomposing 
the power spectrum into two parts as 


Pf{k)=PfAk)+PfMk), (13) 


the hrst and second terms representing the contribution 
from the four- and two-point correlation parts, respec¬ 
tively. Since is related to the power spectrum of 
density fluctuations, P{k), through the Fourier transfor¬ 
mation, the second term, Pf^ 2 (k), is given by 

PfAk) = 2 1 -0^ P{p)Pi\k-p\). (14) 


Note that P{k) is the non-linear power spectrum of den¬ 
sity fluctuations. We calculate it using a halo model 
approach (see also Appendix mi. Figure El shows 
the two-point correlation contribution to the “dimen¬ 
sionless power spectrum” of /, evaluated at z = 0 for 
Mmin = 10 ®M©; the quantity A'j{k) is defined by 


A}{k) 


fc® Pfik) 

2^2 (^2)2 ■ 


(15) 


The “1-halo term” represents correlations between parti¬ 
cles within the same halo, whereas the “ 2 -halo term” rep¬ 
resents correlations between particles in two distinct ha¬ 
los. In other words, two points considered, x and y, are in 
one same halo for the 1 -halo term, and are in two distinct 
halos for the 2-halo term.® The corresponding power 
spectrum will be denoted as P/_ 4 (fc) = Pj\{k) -|- P^'Kk). 
According to Ref. i the 1-halo and 2-halo terms of the 
two-point correlation function are given respectively by 


r(i) 


r(2) 


dn / M 


(17) 


dM ^ ) / dxi 

d Mmin dM \^mPcJ J 

X u‘^{x — Xi\M)u'^{y — Xi\M), 

r dn dn fMiM2 

JMmin ' dM, J^min ^ dM2 V 

.Jdx,Jdx2uH^-x,\m,) 

X u^{y - X2\m2)^l^,l{xi,X2\Mi,M2), (18) 


/o') 

where ^f^fl{xi,X 2 \Mi, M 2 ) is the two-point correlation 
function of halos of mass M, and M 2 , which is given 
by dl^(’’l^i>^ 2 ) « 6 (Mi) 6 (M 2 ) 4 n (’’) HM) is the lin- 

/o') 

ear halo bias and is the linear correlation function]. 
Equations 03 and cn are then Fourier transformed, 
giving the power spectrum as follows: 



X P^“(fc), 



M 

^mPc 


v^{k\M), (19) 

1 2 

b{M)v{k\M) 

( 20 ) 


where v{k\M) is the Fourier transform of 
u'^{x\M)M/ilmPc, which is analytically solvable for 
the NFW profile with fairly a lengthy expression, and 
shown in Fig. El for various values of M. Note that these 
expressions are identical to those for the ordinary matter 
power spectrum, Eqs. m and if we replace 

v{k\M) by u{k\M). 

FigureElshows the dimensionless power spectrum from 
the four-point contribution, A 2 ^(fc), evaluated at 2 : = 0 
for Mmin = 1 O®M 0 . We find that the four-point term al¬ 
ways dominates over the two-point term at all wavenum¬ 
bers; thus, the two-point contribution may be safely ig¬ 
nored. The 1-halo term dominates at smaller spatial 
scales, k > 2 h Mpc“^, whereas the 2-halo term dom¬ 
inates at larger scales. 


2. Four-point contribution 

The four-point correlation function in Eq. m is fur¬ 
ther divided into 1 -halo and 2 -halo contributions, 

i^^\x,x,y,y) +/nL- (16) 


The full expression of con tain s one 1-halo, six 2-halo, seven 
3-halo, and one 4-halo terms 113. However, most of them are 
vanishingly small for a particular configuration considered here, 
03, y, y), as halos are spatially exclusive — the same reason 
as that stated in the second footnote. 
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k [h Mpc 

FIG. 3: The function that determines the four-point contribu¬ 
tion to the power spectrum of CGB anisotropy, v{k\M) [see 
Eqs. 11191 and ISOHi as a function of k evaluated at 2 = 0 
for various values of M. The value of M labeling each curve 
is spaced by one order of magnitude, with the largest and 
smallest values as indicated. 


In Fig. [3 we plot Aj(fc) at z = 0 for the smaller min¬ 
imum mass, The contributions from 

halos whose masses are in some specific ranges are also 
shown. We find an interesting trade-off between large and 
small mass halos. Larger mass halos give much larger 1- 
halo term contributions and thus completely dominate 
Aj(fc) at small scales. On the other hand, smaller mass 
halos give much larger 2-halo term contributions and thus 
completely dominate Aj(fc) at large scales. This phe¬ 
nomena results in a “break” in the shape of Aj(fc) at a 
critical wavenumber, k ^ 1-10 h Mpc~^ depending on 
Afmin, below which the dominant contribution to A^(fc) 
comes from the smallest halos, and above which the dom¬ 
inant contribution comes from the largest halos. 

We find by comparing Figs. 0] and 0 that the 2-halo 
contribution to Aj(fc) is virtually unaffected by Mmin. 
This is somewhat surprising, given that the 2-halo con¬ 
tribution is dominated by the smallest halos having 
Tfmin. The reason for this phenomenon is because we 
define A^(fc) as a ratio of Pf(k) and the clumping fac¬ 
tor squared, (<5^)^ [see Eq. (IT^ j. Reducing the mini¬ 
mum mass increases both Pf{k) and (5^)^ by the same 
amount, which makes the 2-halo contribution nearly in¬ 
dependent of the choice of M^i^. On the other hand, the 
1-halo contribution decreases as Mmin decreases. This 
is easy to understand. As the 1-halo contribution to 
Pf{k) is entirely dominated by the largest halos, it is 



k [h Mpc“'] 


FIG. 4: Dimensionless power spectrum of / = <5^ — (d^) from 
the four-point contribution, A‘j ,^{k) [Eq. lO i. by the 1-halo 
(dotted) and 2-halo (dashed) terms [Eqs. Ill 91 and 121)1 1. Both 
are evaluated at 2 = 0 for Mmin = 10® M©. 



k [h Mpc“'] 


FIG. 5: The same as Fig. |1] but for the smaller minimum 
mass, Mmin = 1O~®M0. Contributions from halos of various 
mass ranges are also shown. 















independent of MminJ however, as the clumping factor 
increases for smaller -Mmin, the dimensionless power spec¬ 
trum, A‘j{k) cx Pf{k)/{6‘^y, decreases for smaller Mmin- 
As a result, a smaller Mj^in makes the shape of A‘j:{k) 
flatter, reducing the power at larger k. A “shoulder” 
due to the 2 -halo term then becomes more prominent for 
smaller Mmin, which is a characteristic that may be ob¬ 
servable in the angular anisotropy, as we discuss in the 
next section. 


IV. ANGULAR POWER SPECTRUM 

By using Eq. m with the 3D power spectrum shown 
in Figs. 0and[31 we can calculate the angular power spec¬ 
trum, Cl, as a function of multipoles, Z, for a given dark 
matter model. We here note that Ci roughly corresponds 
to the correlation between two points on the sky sepa¬ 
rated by an angle 0 « tt/L 

A. Supersymmetric neutralino 

First, we consider supersymmetric neutralinos as the 
dark matter candidate, and assume that the neutralino 
mass is 100 GeV. While the mean intensity, {I-y), is very 
sensitive to the dark matter mass, {I^) oc this fac¬ 

tor exactly cancels out when the amplitude of anisotropy 
is divided by the mean intensity. Thus, Ci depends on 
ruj;. only weakly through the mass dependence of the 
gamma-ray spectrum per annihilation. We evaluate the 
anisotropy power per logarithmic range of I, 

In Figs. ISa) and m, we show the predicted angu¬ 
lar power spectrum evaluated at the observed gamma- 
ray energy of = 10 GeV for Mmin = IO^Mq and 
10~^Mq, respectively. One can understand the shape 
of the angular spectrum for both cases using the same 
argument for understanding the shape of the 3D dimen¬ 
sionless power spectrum, Aj(fc), in the previous section. 
A rapid increase in Ci at small angular scales (large 1) 
is due to the 1-halo term contribution. The 1-halo term 
dominates even at relatively small I, I ~ 10 , for the larger 
minimum mass, whereas it is suppressed significantly for 
the smaller minimum mass, in agreement with the depen¬ 
dence of the shape of Aj(fc) on Mmin that we discussed 
in the previous section. On the other hand, the 2-halo 
term contribution is nearly independent of Mmin, which 
is also in agreement with Aj(fc). (The 1-halo term is 
roughly proportional to the inverse of the clumping fac¬ 
tor squared, while the 2 -halo term is nearly independent 
of the clumping factor.) This peculiar dependence of the 
shape of Ci on opens up an exciting possibility that 
one can “measure” the minimum mass from the shape of 
Cl. Together with the information from the energy spec¬ 
trum of the GGB, therefore, it may be possible to identify 
the mechanism by which gamma rays are produced from 
annihilation of dark matter and the degree to which small 
halos are tidally disrupted by the structure formation. 



1 10 102 103 


Multipole 1 

FIG. 6: Angular power spectrum of the CGB, Ci, from an¬ 
nihilation of supersymmetric neutralinos, evaluated for (a) 
Afmin = 10 ®Mq and (b) Mmin = 1O“®M0. Note that Ci is 
dimensionless: the mean intensity squared should be multi¬ 
plied in order to convert it to the units of intensity squared. 
The neutralino mass rrix is assumed to be 100 GeV. The pre¬ 
dicted angular spectrum is shown at the observed gamma-ray 
energy of — 10 GeV. Gontributions to Ci from the 1-halo 
(dotted) and 2-halo (dashed) terms are shown as well as the 
total signal (solid). 


Is the predicted angular power spectrum detectable? 
We compare the predicted power spectrum with the ex¬ 
pected sensitivity of the GLAST experiment. We take 
the following specifications for GLAST: the field of view 
is Dfov = 47r/fov = 2.4 sr, the angular resolution is 
(Th = 0.115°, and the effective area is Aes = 10^ cm^ 
at 10 GeV |^. Note that the angular resolution is de¬ 
fined as a half width of a Gaussian point spread function, 
and the full width at half maximum is given by •\/81n2crb. 
In addition, for the diffuse gamma-ray observation, the 
background contamination can be reduced to 6 % of the 
CGB, which is a promising characteristic Therefore, 
a fractional error of Ci, 


l 2{\ECna/WfCiY 

Cl V (2Z + l)A;/fov ’ ^ ^ 

is essentially determined by the Poisson noise of the cos¬ 
mic signal. Here A/ is the bin width of I (we choose 
Al = 0.31 for Mmin = IO^Mq and 0.51 for IQ-^Mq), 
Cat ,7 = ^iov[{NN/Nj)'^/N n + l/A’ 7 ] is the power spec¬ 
trum of photon noise, AAr and Afy are the count number 
of backgrounds and the signal {NN/N.y <C 1 for GLAST), 
respectively, and Wi is the window function of a Gaus- 
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sian point spread function, Wi = exp(—Z^cr^/2). If we as¬ 
sume that the CGB detected in GeV region is dominated 
by gamma rays from dark matter annihilation, we may 
use the observed CGB intensity, = 1.5 x 10“® 

GeV cm“^ s“^ sr“^ QIH, to estimate the signal-to- 
noise for anisotropy. The expected number of photons, 
for GLAST would then be = E^{I^)V,{ov^eet = 
10^(t/l yr) at Ej = 10 GeV, while iVjv is negligible. 

In Fig. [TJ we show the predicted angular power spec¬ 
trum at the observed gamma-ray energies of E^ = 3, 
10, and 20 GeV, assuming Mmin = 10®M©, with the ex¬ 
pected la errors of Ci at E^ = 10 GeV for t = 1 yr of 
observations. We find that the GLAST should be able to 
measure the angular power spectrum of the CGB fairly 
easily for 1 year of observations, if the dark matter par¬ 
ticle is the neutralino with mass around 100 GeV and 
its annihilation dominates the observed CGB in GeV re¬ 
gion. The angular power spectrum for the smaller min¬ 
imum mass, Mmin = 1O“®M0 is shown in Fig. O with 
the expected error bars. We find that anisotropy is still 
easily detectable with the GLAST for 1 year of observa¬ 
tions. Therefore, we conclude that, if dark matter par¬ 
ticles are supersymmetric neutralinos and the observed 
CGB in GeV region is dominated by their annihilation, 
the GLAST should be able to measure the angular power 
spectrum of CGB anisotropy, regardless of the minimum 
mass. 

As prospects for detecting CGB anisotropy in the 
GLAST data are very good, it is tempting to specu¬ 
late that one may actually detect CGB anisotropy in 
the existing EGRET data. The EGRET parameters are 
^^fov = 0.5 sr, tTf, = 0.5°, and Aeff = 7 x 10^ cm^ at 10 
GeV and on-axis direction (Apff decreases by a factor of 
2 at 20° and of 6 at 30°) |^. The background noise has 
been suppressed by an order of magnitude compared with 
the photon counts from the CGB. Although the detector 
sensitivity of the EGRET is not as good as that of the 
GLAST, we expect that the Poisson error (photon noise) 
of the signal to be quite small at I < 400 owing to a longer 
integration time (EGRET has observed the sky for much 
more than 1 year). The data above I « 400 would not 
be very useful because of the limited angular resolution. 
While one needs to carefully characterize the systematic 
uncertainty due to the Galactic foreground removal, the 
power spectrum of CGB anisotropy in the EGRET data 
should be measured in search of the signatures of dark 
matter annihilation. 


B. MeV dark matter 

As for the MeV dark matter, the continuum gamma 
rays are emitted by the internal bremsstrahlung, and the 
spectrum shape is shown in Sec. E] (see also. Refs. 

13^1. Figure El shows the predicted angular power spec¬ 
trum of CGB anisotropy for the MeV dark matter with 
= 20 MeV. We show Ci at the observed gamma- 
ray energy of = 3 and 10 MeV, as well as the ex- 



Multipole 1 


FIG. 7: Angular power spectrum of the CGB, Ci, from anni¬ 
hilation of supersymmetric neutralinos, evaluated for Mmin = 
10® Mq. The neutralino mass is assumed to be 100 GeV. 
The predicted angular spectrum is shown at the observed 
gamma-ray energy of = 3, 10, and 20 GeV. The Icr er¬ 
ror bars of Ci expected from GLAST for 1 year of operation 
are also shown at Ej = 10 GeV. 
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FIG. 8 : The same as Fig. Q but evaluated for the smaller 
minimum mass, Mmin = 1 O~®M 0 , which is closer to the free- 
streaming mass of neutralinos. (Note that the precise value of 
the free-streaming mass depends strongly on the dark matter 
mass and interaction properties.) 
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pected uncertainty at 10 MeV. The shape as well as the 
normalization are very similar to those for the neutrali- 
nos. In the MeV regime, the most promising technique 
to detect gamma rays is the Compton scattering inside 
the detector. Although the design of future detectors 
of this kind, such as the Advanced Compton Telescope 
(ACT) [^, are still in progress, we parameterize it with 
tJb = 1 °, Afi ff = 100 cm^, Ofov = 1 sr, t = 1 yr, and 
AOv/Aly = 1 [^. Again taking the mean CGB intensity 
at 10 MeV to be what has been observed, 10“"^ cm“^ s“^ 
sr“^ [a Elia, the signal count with ACT would be 
Nj = 3 X 10®(t/l yr). With these parameters for ACT, 
we find that it should be easy to measure the power spec¬ 
trum in the range of 10 < ^ < 100. Above I = 100, the 
errors become exponentially larger because of the detec¬ 
tor angular resolution, i.e., due to the effect of the window 
function Wi in Eq. 121); this is unavoidable unless better 
angular resolution is realized. 

The dependence of Ci on the minimum mass, Mmin, is 
almost the same as that for the neutralino. In Fig. nni we 
show the angular power spectrum for the smaller mini¬ 
mum mass, Minin = 1 O“®M 0 . (And this is a correct 
mass to use, up to the dependence on the dark mat¬ 
ter mass 1 ^, when the emission is due to the inter¬ 
nal bremsstrahlung; for the minimum mass appropriate 
for the 511 keV emission, see arguments in Ref. 
Again, the predicted angular power spectrum is well 
above the expected uncertainty, and thus our conclusion 
about prospects for detection of anisotropy is robust re¬ 
gardless of the minimum mass. 

The COMPTEL has measured the CGB in this energy 
region [3- The parameters of the COMPTEL are flfov = 
1 sr, ab = 1.5°, Aeff = 30 cm^ at 20 MeV which are 
almost comparable to the parameters adopted for the 
ACT. The biggest issue is, however, its poor signal-to- 
background ratio, Nj/Nn « 0 . 01 ; thus, it would be quite 
challenging to detect CGB anisotropy in the COMPTEL 
data because of large statistical errors. 


V. DISCUSSION 

A. Dependence on gamma-ray energy 

In addition to the shape of the angular power spectrum 
as a function of Z, it would be useful to discuss dependence 
on the gamma-ray energy in order to obtain further infor¬ 
mation on the CGB origin. We show the energy spectrum 
ofisotropic, {Ij), and anisotropic, {Ij)[l{l+l)Ci/{2 t:)Y/^ , 
components of the GGB in units of intensity in Fig. da), 
for both the neutralino (m^ = 100 GeV) and MeV dark 
matter (m^ = 20 MeV) assuming Mmin = 10 ®Mq. The 
energy spectrum of anisotropy divided by the mean in¬ 
tensity, l{l + 1 )Ci/{2tt), is also shown in Fig. ITTT bL We 
show only I = 100. The isotropic spectrum is boosted ac¬ 
cordingly in order to explain the observed intensity; the 
details are given in the last paragraph of Sec. 11. 

We find that anisotropy is enhanced when the spec- 



FIG. 9: The same as Fig. Q but for MeV dark matter with 
nix = 20 MeV, evaluated at = 3 (dotted) and 10 (solid) 
MeV. Error bars are calculated for ACT in 1 year of operation. 
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FIG. 10: The same as Fig. El but evaluated for the smaller 
minimum mass, Mmin = 1 O~®M 0 , which is closer to the free- 
streaming mass of MeV dark matter. (Note that the precise 
value of the free-streaming mass depends strongly on the dark 
matter mass and interaction properties.) 
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E, [GeV] 

FIG. 11: (a) Energy spectrum of isotropic (dashed) and 

anisotropic (solid) component of the CGB intensity for = 
100 GeV (right curves) and 20 MeV (left curves), evaluated 
for Mmin = 10 ® Mq. Anisotropic components are evaluated 
at I — 100. (b) Energy spectrum of anisotropy divided by the 
mean intensity at I = 100. 


trum has sharp features such as edges (as seen at ^ 20 
MeV for the MeV dark matter case) or lines. This has 
been pointed out by Ref. in the context of CGB 
anisotropy from Type la supernovae. Therefore, if the 
energy resolution of gamma-ray detectors is sufficiently 
good to resolve these properties, the spectral features 
would provide another powerful diagnosis of dark mat¬ 
ter annihilation signals. The features might also al¬ 
low us to obtain information on the redshift evolution 
of the sources (see. Ref. for more details). Other 
potential spectrum signatures include the line gamma 
emission (xx XL ^*^7) internal bremsstrahlung 
ixx W~W^^ 1^) from annihilation of neutralinos, 
as well as the 511 keV line emission from annihilation of 
MeV dark matter particles into pairs. 


B. Other astrophysical sources 

The angular power spectrum shown in Figs. [71177)1 
should be very characteristic of annihilating dark matter 
in extragalactic dark matter halos, as the gamma-ray in¬ 
tensity is proportional to the density squared. The inten¬ 
sity of gamma rays coming directly from other astrophys¬ 
ical sources should be linearly proportional to density. It 
is likely that blazars are the most dominant constituent 
of the GeV gamma rays over a wide energy range [zaEii 
Assuming that blazars are biased tracers of the underly¬ 


ing mass distribution, the two-point correlation function 
of blazars should be simply given by that of density fluc¬ 
tuations, P{k). For more quantitative study, we perform 
the following simple analyses for the power spectrum of 
blazars. 

First, we assume that blazars are very rare objects and 
their angular spectrum is entirely dominated by the shot 
noise. In this case the angular power spectrum does not 
depend on I, and thus l{l + l)Ci is proportional to P at 
I ^ 1. We show in Figs. lEKa) andlEKb) that the shot 
noise spectrum totally lacks the power on large angular 
scales (i.e., the spectrum is too steep), and can be easily 
distinguished from dark matter annihilation. 

Second, we take the other extreme limit where blazars 
are quite common and trace the underlying matter den¬ 
sity field, 6, fairly well. Specifically, we calculate the 
line-of-sight integral in Eq. m with S'^ replaced by 6. 
The volume emissivity is then given by P.y oc i5(l-|- 
z)^EjdN^/dEj, and we assume a power-law blazar en¬ 
ergy spectrum falling off as dN^/dE^ oc E~'^ (e.g., 
for Mrk 421), which gives W{E^,z) oc We 

calculate the angular power spectrum, Ci , using Eq. m 
with Pf{k) replaced by P{k). Finally, we assume that 
blazars form only when the mass of host dark matter 
halos is larger than Mmin = 10 ^^Mq. In this simplified 
prescription, the average number of blazars in a halo lin¬ 
early increases with the host halo mass. In reality, how¬ 
ever, this may not be true and we may also need to take 
into account the difference between a central galaxy and 
satellite galaxies within a dark matter halo. This kind 
of galaxy distribution model, called the Halo Occupation 
Distribution (see Ref. m for a review), could be made 
more realistic for blazars; however, we do not pursue it in 
any more detail, and adopt the above model for calculat¬ 
ing the angular power spectrum of gamma-ray emission 
from “blazars.”^ 

We show the shape of the angular power spectrum ex¬ 
pected from blazars in Figs. ll2r a) and ll2r b). (Normaliza¬ 
tions are taken arbitrarily.) The blazar power spectrum 
is determined by the 1-halo term at all multipoles; thus, 
the spectrum is nearly a power law with a little bit of flat¬ 
tening above I ^ 200, which is due to the minimum mass 
of halos hosting blazars, 10^^Mq. The small-scale power 
drops as we raise the minimum mass of host halos (see, 
e.g.. Fig. 2 of Ref. [i^). In Fig. IT^ a). we compare the 
blazar spectrum with the spectrum of GGB anisotropy 
from annihilation of neutralinos with = 100 GeV for 
the large minimum mass, 10®Mq. As we have noted be¬ 
fore, the annihilation spectrum for M^nin = 10®Mq is also 


^ Our model is probably not a very good model for describing the 
power spectrum of real blazars. It simply describes the angular 
power spectrum of dark matter halos above which emit 

gamma rays (by whatever mechanism) with an energy spectrum 
falling off as dNj/dEn/ oc . In reality, not all dark matter 
halos above 10 ^^Mq host blazars, and blazars are known to be 
transient objects which turn on and off. 
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FIG. 12: Shape of the angular power spectrum of the CGB 
expected from unresolved “blazars” (solid lines; see the foot¬ 
note 5 for the reason why the quotation marks are put) with 
arbitrary normalizations. The power spectrum from anni¬ 
hilation of neutralinos with = 100 GeV is also plot¬ 
ted as the dashed lines. The adopted gamma-ray energy is 
10 GeV, and the minimum mass of dark matter halo is (a) 
10®Mq, and (b) 10“®Mq. The dotted lines show the shot 
noise {C'l = const.) with arbitrary normalizations, which rep¬ 
resent the power spectrum of very rare sources. 


dominated by the 1-halo term at all multipoles, which 
makes the annihilation and blazar spectrum actually look 
similar at I < 200. Above I ~ 200, however, the annihila¬ 
tion spectrum continues to grow whereas the blazar spec¬ 
trum flattens out. In Fig. mb), we use the smaller min¬ 
imum mass, which is probably more realistic 

for dark matter annihilation. In this case the dominant 
contribution comes from the 2-halo term, which makes 
the annihilation spectrum substantially flatter than the 
blazar spectrum. In other words, the annihilation spec¬ 
trum has much more power at large angular scales, which 
should be easily distinguished from the blazar spectrum. 
The same argument should also apply to the CGB in the 
MeV region from Type la supernovae. (See Ref. for 
another approach to calculating the angular correlation 
function of Type la supernovae.) 

Another potential astrophysical contributor to the 
CGB is the gamma-ray emission from particles acceler- 
ate d by shock waves forming in galaxy clusters III El 
El El El El’ There have been a few claims about 
detection of a correlation between the position of clus¬ 
ters and the EGRET data IH ■ When the gamma 

rays are emitted mainly by the inverse Compton scat¬ 
tering off the CMB photons by accelerated electrons. 


which many authors considered to be the dominant emis¬ 
sion process, the gamma-ray intensity would be propor¬ 
tional to the baryon density at the shock. On the other 
hand, when the gamma rays are produced by collisions 
between cosmic-ray protons (accelerated by the cluster 
shocks) and intracluster medium, the gamma-ray inten¬ 
sity would be proportional to the baryon density squared 
at the shock. However, the baryon density distribution 
at the shock would be very different from the distribution 
of dark matter particles. One would need cosmological 
hydrodynamical simulations to study the shape of the 
angular power spectrum of CGB from such processes. 


C. Substructure of dark matter halos 

The biggest uncertainty in calculations of the CGB, 
whether the mean intensity or anisotropy, is the substruc¬ 
ture (i.e., small-scale dumpiness) within dark matter ha¬ 
los. A^-body simulations of collisionless dark matter uni¬ 
versally indicate the existence of sub-halo clumping of 
dark matter . This effect has been entirely ignored in 
our calculations, as we have assumed that each halo has 
a smooth density profile. It is expected that the dark 
matter substructure increases the gamma-ray intensity 
by simply increasing the dumpiness.® How much it is in¬ 
creased depends very much on the exact amount of sub¬ 
structure, which is difficult to predict due to a limited 
resolution of the current A^-body simulations. In this 
regard, our calculations provide the lower limit to the 
angular power spectrum of CGB. The fact that our pre¬ 
dictions are already well above the noise level of GLAST 
and ACT is extremely encouraging for prospects of de¬ 
tection of anisotropy. 

Let us comment on some effects of dark matter sub¬ 
structures on the angular power spectrum. It is ex¬ 
pected that the effect of substructure would be smaller 
for anisotropy than for the mean intensity, when the am¬ 
plitude of anisotropy is divided by the mean intensity 
(which is what we have been using to calculate Ci). But 
still, the increased clumping will increase the four-point 
function of density fluctuations, particularly 1-halo term 
which dominates at small scales. Therefore, we expect 
the substructure to modify Ci at large I, while we do 
not expect a significant change in Ci at small I, where 
the four-point function is dominated by the 2-halo term. 
Finally, in the limit that dark matter annihilation sig¬ 
nals are completely dominated by the sub-halos within 
a bigger (host) halo, e.g.. Earth-mass micro halos in a 
Milky-Way size halo, the gamma-ray intensity may be 
linearly proportional to the density profile of host halos 
|2f)| . In this model, there is no gamma-ray emission from 
the intra-halo medium between the sub-halos, and the 


For gamma rays from annihilatio n of dark matter in the sub- 
structure of our Galaxy, see Refs. 1^ IsS. 1^ IsS. 1^ IsSi . 
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sub-halos are assumed to follow the density profile of a 
host halo. Our formalism should be modified for this case 
by including the number of sub-halos within a host halo 
as a function of the host halo mass, which is essentially 
the Halo Occupation Distribution for the sub-halos. 

Incorporating dark matter substructures into the halo 
model is an active field of research, and therefore we ex¬ 
pect to be able to get a better handle on this effect in 
near future. 


VI. CONCLUSIONS 

The CGB provides indirect means to probe the nature 
of dark matter particles via high-energy photons from 
dark matter annihilation. The dark matter annihilation 
has occurred in all the past halos, and now contributes 
to the CGB flux at some level; its contribution might be 
dominant if the flux is boosted by some mechanism such 
as tidally survived dark matter substructure. Therefore, 
revealing the origin of the CGB is an very important 
problem that is potentially connected to the dark mat¬ 
ter properties and halo substructure as well as ordinary 
astrophysical objects. In order to achieve this purpose, 
while the energy spectrum of the mean intensity of CGB 
has been investigated by many researchers, anisotropy 
of CGB from dark matter annihilation has been entirely 
neglected. In this paper we have calculated the angu¬ 
lar power spectrum of CGB from dark matter annihila¬ 
tion for the first time, using an analytical halo model 
approach. As for dark matter candidates, we have dis¬ 
cussed two possibilities: one is the supersymmetric neu- 
tralino, one of the most popular candidates today, which 
contributes to the CGB in GeV region. The other is the 
MeV dark matter, a dark matter species first introduced 
to explain the 511 keV emission line from the Galactic 
center by annihilation into e“e“'" pair, which contributes 
to the CGB in MeV region. 

Since the gamma-ray intensity from annihilation is pro¬ 
portional to the density squared, we have derived 
the 3D power spectrum of the quantity / = — (^^), 

Pf{k). The calculation involves the Fourier transforma¬ 
tion of the four-point correlation function of underlying 
mass density fluctuations, which has been shown 

to dominate over the two-point correlation contribution. 
This Pf_i(k) includes the 1-halo and 2-halo terms, 
the former containing two points within the same halo, 
and the latter containing two points in two different ha¬ 
los. The analytical expressions for Pf{k) from each term 
are given in Eqs. ra . ira . and CT . and the results are 
plotted in Figs. [3 El and El At all scales the four-point 
contribution totally dominates the signal; at small scales, 
the 1-halo term dominates. 

We note that our formalism can also be used for any 
other emission processes that involve collisions of two 
particles. For example, one may use this to compute 
the power spectrum of free-free or bound-free emission 
from the ionized gas in halos. For this application one 


needs to replace the dark matter density profile with the 
gas density profile. 

Using Eq. m that connects the 3D power spectrum, 
Pf{k), to the angular power spectrum, Ci, we have cal¬ 
culated the GGB angular power spectrum as a function 
of multipoles, I, for both the neutralino = 100 GeV) 
and MeV dark matter (m^ = 20 MeV) at various gamma- 
ray energies. We have also compared the predicted sig¬ 
nals with the expected sensitivity of future gamma-ray 
detectors — GLAST in GeV region (for neutralinos) and 
AGT in MeV region (for MeV dark matter). The re¬ 
sults are shown in Figs. CHini For both cases, we have 
found that these detectors will have sufficient sensitivity 
to measure the angular power spectrum with reasonable 
accuracy. 

We have studied the effects of the minimum mass, 
Mmin, on the predicted angular power spectrum in de¬ 
tail. While the 1-halo contribution, which dominates at 
small angular scales (large 1), decreases for smaller Mmin, 
the 2-halo contribution, which dominates at large angu¬ 
lar scales (small 1), is virtually unaffected by Mmin. This 
property results in a peculiar dependence of the shape 
of Cl on Mminj which may be used in combination with 
the information from the energy spectrum of the GGB to 
determine As Mmin depends on the radiation pro¬ 

cesses of gamma rays from annihilation as well as the sur¬ 
vival of micro halos contributing to the GGB, the shape 
of Cl provides a powerful tool for determining these prop¬ 
erties, which are otherwise difficult to probe. Our conclu¬ 
sion about prospects for measuring Ci of GGB anisotropy 
are robust regardless of M^i^- 

By applying our formalism with some modification, 
we have shown that the other astrophysical sources such 
as blazars would reveal a different shape of the angu¬ 
lar power spectrum, as these contributions are linearly 
proportional to density fluctuations. The shape of Ci at 
large I might further change when we take into account 
the dark matter substructure, but the result in this pa¬ 
per sets a lower bound on anisotropy at all multipoles, 
which provides excellent prospects for detection of GGB 
anisotropy by future gamma-ray detectors. 

We conclude that the angular power spectrum of GGB 
provides a smoking-gun signature for gamma-ray emis¬ 
sion from annihilation of dark matter particles, which 
would be a powerful tool for understanding the nature of 
dark matter particles. 
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APPENDIX A: RELATION BETWEEN 
ANGULAR AND THREE-DIMENSIONAL 
POWER SPECTRUM 


In this section, we derive the relation between the an¬ 
gular power spectrum, Ci, and the 3D power spectrum, 
Pf{k), that we have given in Eq. (ITnil . With the spatial 
Fourier transformation, Eq. is rewritten as 

= J dnj dr J fk{r)p'"'^W{r)Yi*^{n) 

^ I I I 

X 31'{kr)Y{:^,{k)Yi,rn'{fi) 

I'm' 

= j drW{r) J ^ fk{r)ji{kr)Yi*^{k), 

(Al) 

where in the second equality is expanded into the 
spherical harmonics and the spherical Bessel function, 
ji(kr), using Rayleigh’s formula, and the last equality is 
due to the orthonormal relation of Yim{n). Hence, the 
angular power spectrum Ci = (|aimP) is expressed by 


{I^fCi = J dr W{r) J dr' W{r')^ J dk Pf{k;r,r') 

X ji{kr)ji{kr')Yim{k)Yir^,{k) 

= J dr W{r) J dr' W(r') 


k^dk Pf{k;r,r')ji(kr)ji{kr') 


(A2) 


where we used the definition of P/, (/fc(r)/^,(r')) = 
{2TT)^S^^\k — k')Pf{k]r,r'). As a final step, using the 


following approximation, 

— J k^dk Pf{k-,r,r')ji{kr)ji{kr') 

which is valid if Pf{k',r,r') varies relatively slowly as a 
function of /c, we arrive at Eq. ca. 


APPENDIX B: POWER SPECTRUM OF 
DENSITY FLUCTUATION 


In this section, we summarize the power spectrum P{k) 
of density fluctuations, S. The detailed derivation has 
been given in, e.g., Refs. minEi , to which we refer 
the reader. P{k) includes two distinctive terms, P{k) = 
P^^{k) -I- P‘^^{k), each term given by 




Wr. 


\u{k\M)\^ 

b{M)u{k\M) 


(Bl) 
1 2 


X Plin(^); 


(B2) 


representing the 1-halo and 2-halo contributions, respec¬ 
tively. Here u{k\M) is the Fourier transform of u{r\M), 
the analytic representation for the NEW profile given by 
Eq. (81) of Ref. is the bias parameter and 

Piin is the linear power spectrum, for which we adopt the 
fitting formula in Ref. [^. The values for the 1-halo 
and 2-halo terms as well as the linear power spectrum 
are shown in, e.g., Fig. 1 of Ref. |^. The 1-halo term 
dominates at small spatial scales (large k). 
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